Full-counting statistics for molecular junctions, 
the fluctuation theorem and singularities 



Y. Utsumi, 1 ^ 0. Entin-Wohlman, 2 '[t| A. Ueda, 3 and A. Aharony 2 ^ 

1 Department of Physics Engineering, Faculty of Engineering, Mie University, Tsu, Mie, 514-8507, Japan 
2 Department of Physics and the Use Katz Center for Meso- and Nano-Scale Science and Technology, 
Ben Gurion University, Beer Sheva 84105, Israel 
1 Faculty of Pure and Applied Sciences, Division of Applied Physics, 
University of Tsukuba, Tsukuba, Ibaraki, 305-8573, Japan 
(Dated: October 9, 2012) 

We study the full-counting statistics of charges transmitted through a single-level quantum dot 
weakly coupled to a local Einstein phonon which causes fluctuations in the dot energy. An analytic 
expression for the cumulant-generating function, accurate up to second order in the electron-phonon 
coupling and valid for finite voltages and temperatures, is obtained in the extended wide-band 
limit. The result accounts for nonequilibrium phonon distributions induced by the source-drain 
bias voltage, and concomitantly satisfies the fluctuation theorem. Extending the counting field to 
the complex plane, we investigate the locations of possible singularities of the cumulant-generating 
function, and exploit them to identify regimes in which the electron transfer is affected differently 
by the coupling to the phonons. Within a large-deviation analysis, we find a kink in the probability 
distribution, analogous to a first-order phase transition in thermodynamics, which would be a unique 
hallmark of the electron-phonon correlations. 
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I. INTRODUCTION 

It has been revealed in electric transport measurements 
that the electron-phonon interaction induces unique fea- 
tures in the nonequilibrium current through molecular 
junctions^ and atomic wires^— In particular, when 
the phonon energy is small compared to the resonance 
width on the junction, the inelastic phonon scattering 
increases/decreases the current for small/large transmis- 
sion probabilities as the source-drain bias voltage V ex- 
ceeds the local phonon frequency ujq^-~— Quite a number 
of theoretical microscopic models, e.g. Ref. 0, including 
density-functional theories; 7 - have been devoted to such 
junctions (see Ref. § and references therein). They re- 
vealed that a simplified model^ of a single-level quantum 
dot coupled to a local Einstein phonon mode causing fluc- 
tuations of the dot energy level^ seems to suffice to cap- 
ture this ubiquitous feature. Recently, the current noise 
of an atomic wire has been measured^ and it was ob- 
served that the electron-phonon interaction can enhance 
or reduce the noise, depending on the value of the trans- 
mission probability. Based on the theory of Avriller and 
Levy Yeyatijii the negative correction is understood as 
resulting from the anti-bunching of two electrons: An 
electron cannot be inelastically scattered by a phonon 
when the final state is already occupied by another elec- 
tron. 

Avriller and Levy Yeyati considerations follow from 
the theory of full-counting statistics (FCS)) 12 i 13 which 
is most convenient for analyzing nonequilibrium electric 
transport. Indeed, considerable effort has been invested 
in recent years in exploiting FCS to study various aspects 
of nonequilibrium quantum transport [e.g. Refs. Illl424l 
and references therein]. FCS refers to the probability 



distribution P T (g) of the charge q transmitted through a 
quantum conductor during a certain measurement time 
t at out-of-equilibrium conditions (we set e = h = 1). 
The effect of coupling to a vibrational mode on the elec- 
tric transport has also been investigated in this context, 
both for a wea k 11 ' 19 " — and a stronger— electron-phonon 
coupling. 

In general, it is a rather formidable task to calculate 
the FCS of interacting electrons. For this reason, most 
of the investigators have taken advantage of the Kcldysh 
field-theory technique. There, the characteristic function 
(CF), 



Z r (A) = ^P T ( 9 )e 



iqX 



(i) 



or the scaled cumulant-generating function (CGF)2£ per- 
taining to the steady state, 

J-(A) = lim -inZJX) , (2) 

r — >oo T 

can be formally written as a 'partition function' or a 'free 
energy', respectively, defined along the Kcldysh time- 
contour. The A appearing in these formulae is termed the 
counting field or the counting parameter. However, even 
if one calculates the CGF successfully, one still needs to 
find a way to characterize the electronic correlations in it. 
One promising approach would be to utilize the location 
distribution of the zeros of the CF, or equivalently that 
of the singularities of the CGF, by allowing the counting 
field A in Eq. ([I} to attain complex values,^— similarly 
to the Yang-Lee theory of phase transition^ This idea is 
based on the recent observation that, upon transforming 
A into u. 



(3) 
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the singularities of the CGF of nonintcracting electrons 
transported between two terminals are all on the neg- 
ative real axis of the u— plane . 26 ' 31 ' 32 It suggests that 
singularities off the negative real axis would character- 
ize electronic correlations. From this aspect, molecular 
junctions are rather advantageous since second-order per- 
turbation theory in the electron-phonon coupling would 
capture most relevant features of the electron-phonon 
correlations in them, allowing for obtaining the location 
distribution of the singularities. 

Another recent ingredient is the fluctuation theorem 
(FT) £2r— The FT is a consequence of micro-reversibility 
and can be understood as a microscopic extension of the 
second law of thermodynamics. Despite its simple ap- 
pearance, a detailed-balance like relation [see Eq. (|6]) 
below], the FT reproduces the linear- response results, 
i.e. it ensures the fluctuation-dissipation theorem and 
Onsagcr's reciprocal relations close to equilibrium^ - — 
while conveying invaluable information at noncquilib- 
rium conditions. For molecular junctions, the FT has 
been addressed using the master-equation approach of 
FCS for incoherent electron transport] 22 ' 23 The FT is 
considered to be a basic symmetry, such as gauge in- 
variance, which the CF should fulfill. 

In the present paper we investigate the FCS of elec- 
trons coupled to phonons under out-of-equilibrium con- 
ditions. Our calculation is accurate up to second order in 
the electron-phonon coupling. We begin in Sec. [TTI with 
brief general explanations of the FT, the large-deviation 
analysis and singularities of the CGF. We then introduce 
in Sec. Mil the model Hamiltonian and present analyti- 
cal results for the CGF and detailed explanations of the 
calculations and the approximations involved. A general 
formulation based on the nonequilibrium Luttinger-Ward 
potential^— is relegated to Appendix lAl In Sec. IIVI we 
analyze the singularities of the CGF and demonstrate the 
probability distribution within the large-deviation anal- 
ysis. Section fVl summarizes our results. 



II. THE FLUCTUATION THEOREM AND 
SINGULARITIES 



The definition of the probability distribution P T (q) in 
quantum systems requires special care . 12 ' 13 Full-counting 
statistics theory begins with introducing the CF [see 
Eq. ([TJ and Eq. (fT4]) below for the definition employed 
in the Keldysh technique] and then defining the quasi- 
probability distribution by the inverse Fourier transform 
of the CF, 



Eq. ©, 



(4) 



More details are given in Sec. Mil The n-th cumulant, in 
steady state, is given by the n-th derivative of the CGF, 



an 



lim 

T— >OQ 



«<?")) d n F(\) 



(5) 



For example, the average current is the first cumulant, 
and the current noise is the second one. 

The FT in the context of quantum electric trans- 
port relates the probability of the current to flow oppo- 
sitely to the bias voltage because of thermal agitations, 
P T (-q)^i-ia to the distribution P T {q), 



P T (-q) = P T (q) 



-ftqV 



(6) 



where (3 is the inverse temperature. The FT can be equiv- 
alently written in terms of the CGF J 7 , Eq. @, as 



T(X) = J r (-X + i/3V). 



(7) 



The relation ([7]) restricts the possible locations of the 
singularities of the CGF in the A— plane. As an ex- 
ample, we depict in Fig. [T] (a) the branch cuts corre- 
sponding to continuous singularities of the CGF per- 
taining to two-terminal transport of nonintcracting elec- 
trons, Eq. below. This CGF is 2n— periodic along 
the real axis of A, which guarantees integer values of 
charged The branch cuts, depicted by thick lines, are 
at Re A = (2n — l)ir, where n is an integer. The FT en- 
sures that the branch cuts are symmetrically distributed 
around A = ifiV/2, and their positions on the upper pos- 
itive/negative real and lower negative/positive real sides 
of the dotted line ImA = (3V/2 in Fig.Q](a) are related. 
The 2n— periodicity is removed by the conformal trans- 
formation Eq. Q . Then the branch cuts are on the neg- 
ative real axis of the u-plane [Fig. Q] (b)] 2*Li21£L21 

The steady-state probability distribution, beyond the 
central-limit theorem, is derived within the theory of 
large deviations^ At steady state, realized in the t^co 



limit, we scale q = It and Z T 



Then the integral of 



the inverse Fourier transform Eq. (TjQ) can be estimated by 
the saddle-point approximation^ and the result is writ- 



ten with the rate function^ 1 as P- 



-tI 



Since P T 



is real and positive, the saddle point is expected to re- 
side on the imaginary axis of complex A— plane. Then the 
rate function is written in the from of a Legendrc-Fcnchcl 
transform, 

1(1) = - lim - lnP.(Ir) = m&x{iXI - J"(A)} . (8) 

t— >oo T A 

Here A is a purely imaginary number. In most cases, the 
CGF is real, i.e. the imaginary part of the exponent of 
the Fourier integral (Tj[]), lnZ r (A) — iXq, is zero on the 
imaginary axis of A— plane. Then the imaginary axis is 
expected to be a steepest path of the integral (j4])4^ There 
are few exceptions where singularities of the CGF are on 
the imaginary axi o 15 ' 16 as we will also find later. The 
relation between the CGF and the rate function is anal- 
ogous to that of thermodynamic potentials. It suggests 
that singularities on the imaginary axis of the A— plane 
would also cause characteristic features in the rate func- 
tion. 



3 




FIG. 1: Singularities of the CGF of noninteracting electrons, 
Eq. (|36[1 . in the A— plane (a) and in the u— plane (b). The 
FT relates the positions of the singularities on the upper pos- 
itive/negative real and lower negative/positive real sides of 
the dotted line (a). In it— plane the FT relates the positions 
of the singularities inside and outside the dotted circle (b). 



III. MODEL HAMILTONIAN AND 
CUMULANT-GENERATING FUNCTION 

Our explicit calculations are carried out for a simple 
model [Fig. [5] , a single- level quantum dot coupled to a 
local Einstein phonon, which induces fluctuations in that 
level£ ~ 10 i 49 The model Hamiltonian is 

^ = ^mol + ^ph + ^load + ^tun,± ' ( 9 ) 

where the ± subscript is due to the presence of counting 
fields (see below). The molecular bridge or the atomic 
wire is described by the Hamiltonian H m oi, 

U mol = [e +7(b + tf))clc , (10) 

in which Co (cj) destroys (creates) an electron on the lo- 
calized level representing the molecule, of energy eo, b and 
&t are the destruction and creation operators of the vibra- 
tions to which the electron is coupled while residing on 
the dot, and 7 is the strength of the electron-phonon cou- 
pling. The vibrational modes obey the harmonic Hamil- 
tonian 

U ph = u b% , (11) 

where ujq is the frequency of the Einstein phonon. The 
leads are represented by free electron gases, of creation 
and destruction operators cj fc and c lk , and eigen energies 

^lcad = E E C rfc C L C rfc ■ ( 12 ) 
r=L,R k 

Finally there is the tunneling Hamiltonian, coupling the 
leads to the molecule. This part of the Hamiltonian is 
augmented by the counting fields. Those appear as phase 
factors on the operators c rk and cj fe , 

?W -EE J r e±iAr 4-Co + H.c. , (13) 

r=L,R k 



FIG. 2: Schematic description of the single-level quantum dot 
coupled to an Einstein vibrational mode. 

where Jl and Jr are the tunneling amplitudes between 
the left and the right lead, and the dot. 

In terms of the tunneling Hamiltonian Eq. (|13p . the 
characteristic function isH 

ZtW = (r K CX P (-i J dtHtun,±{t)l)) > ( 14 ) 

where K denotes the Keldysh contour, which runs from 
t = —t/2 to t = t/2 on the upper branch and returns to 
t = —t/2 on the lower one (see Fig. [3]) and Tk is the time- 
ordering operator along that contour. The superscript I 
indicates time dependence in the interaction picture. The 
± notation indicates the branch of the Keldysh contour 
on which the tunneling Hamiltonian is effective, + for 
the upper branch in Fig. [3] and — for the lower one. In 
the long-time limit, the scaled CGF depends only on the 
the difference of the two counting fields 

A = A L - A R . (15) 

As Al+Ar counts the number of electrons flowing into the 
dot, the fact that the CGF depends solely on A implies 
current conservation. 

K+ 



K- 

t = t/2 t = -t/2 

FIG. 3: The Keldysh contour. 

The FT, including current conservation as represented 
by Eq. (fl~5"j). has been proved using the perturbation ex- 
pansion in the interaction^ Although the proof was for 
Coulomb interaction case, it may be extended to the 
electron-phonon interaction case as well. We even ex- 
pect it to be correct for any closed diagram in the linked- 
cluster expansion (for example, the diagrams in Fig. |4]) 
which suggests, in turn, that by simply summing up 
closed diagrams one will obtain an approximate generat- 
ing function satisfying the FT and current conservation. 

To carry out the perturbation expansion in the present 
case, we utilize the noncquilibrium Luttinger-Ward 
potential^ -47 ' 50 <J>. This theory relics on the self energy 



4 



of the Green function and requires self-consistent calcu- 
lations to ensure various conservation laws. For our case, 
since the electron-phonon coupling 7 is weak, we perform 
the calculation up to second-order accuracy C(7 2 ) (see 
Appendix O for details) 

T(X)= J- (A)-$< 2 >(A) , (16) 

where J"o is the CGF pertaining to noninteracting elec- 
trons, and is given in terms of the Keldysh Green function 
gx of the electronic part of the Hamiltonian Eq. ((9|), 

W) = ^1 dwlndet^M] , (17) 

where the trace is performed over the 2x2 Keldysh space. 
The scaled Luttinger-Ward potential, $( 2 ), appearing in 
Eq. (fl6|) is accurate up to 0(7 2 ), 

lim = #(A) = i>^(A) +0( 7 4 ) . (18) 

r— >oo 7" 



Here is a subtle point. A naive second-order perturbation 
theory is not capable of producing the correct noncqui- 
librium phonon distribution,—^ ' 19 ' 51 as has been recog- 
nized before for a local moment coupled with a nonequi- 
librium electron g&s£^ Therefore, we will first perform 
a re-summation of linked diagrams for the cumulant- 
generating function, which is known to suffice for a 
proper description of the steady stated and then cal- 
culate accurately within second order in the electron- 
phonon coupling. 



A. The Electronic part 



The inverse Keldysh Green function g\ reads 



,9a M~ 



i£ r=L) Rr r [/ r -M-/ r + M]/2 

-*Er = L,R r r/rMe- iA - 



i£r=L lR r r /+( W )e^ 

- w + e + 1 £ r=L)R r r [/-( w ) - /+M]/2 



(19) 



where IY and Tr are the partial widths of the localized 
level induced by the coupling with the leads, 



1-kvAJA 



L,R 



(20) 



Here and vn are the densities of states in the left 
and right leads. Each of the leads is specified by its 
chemical potential /i r , such that /il — Mr = an d their 
electron/hole distribution is accordingly given by 



1 



e ±f)(w-n r ) _|_ I 

The Green function itself is then 



3a M 

and its four components are 



ff ++H g+_-(u>) 
9\ + ( UJ ) g x (w) 



(21) 



(22) 



3a» 



1 



1 



fl x (uj) \ s(u - e ) + iT/2 



r=L.R 



fix 



'(«) = E 



r=L,R 



(23) 



where s = — /+ for s = +/— . The lesser and greater 
Green functions, in the absence of the counting field, g* a , 
are expressed with the density of states on the localized 



level normalized by T — 5Z r=L r the width of the 
resonance on the localized level, 



r 2 /4 



( W -eo) 2 +r2/4' 



as 



(24) 



(25) 



The dependence on the counting field is contained in the 
function 

O / \ det 3A( W ) _1 -1 1 ryi \ 

^aM = - , I / x i = 1 + 7( w ) 
det3 (w) 

x [f+(u;)f-(uj)(e iX 1) + /+H/ L -( W )(e- iA - 1)] . 

(26) 

Here appears the frequency-dependent transmission of 
the localized level, 

7» = ap(w) , (27) 
with the transmission probability at resonance, 



4r L r R 
r 2 



(28) 



Our calculation is confined to the extended wide-band 
limit^ in which the level broadening T or the dot level 
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| eo | is larger than all other energy scales including the 
phonon frequency cj , i- e. 



or 



\V\,k B T,uJo <r, 



|F|,fc B T ! cj ,r< I e 1 



(29) 



(30) 



With this condition, the frequency-dependent normalized 
density of states Eq. (|24[) is replaced by its value at the 
Fermi energy, 



p{uj) ps p Q = 



4e 2 



r 2 ' 



(31) 



Within the extended wide-band limit approximation 
the transmission becomes energy independent, and con- 
sequently the computation of the integral determining 
the zeroth-order CGF [see Eq. (|17|) ] is straightforward. 
The key observation is that the variable transformation 
z = cxp[/3(u/ + (pl + Mi?)/2)] transforms Eq. ([26]) into a 
simpler form, 



(z - Z x+ )(z - Z\J) 



where 



Z x± =X X ± slX( - 1 



and 



X\ =(T — 1) cosh — - — Tcosh ^ ^ — — 



Then we obtain 

dT (X) _ -1 2d lX X x 



d(i\) 2n/3 Z M 



hi 



Z\ — 



Z X — z 



(32) 



(33) 



(34) 



(35) 



and consequently the CGF itself is derived^ from Eq. 
(|35p by integrating over A, 



^o(A) 



(arccoshA A ) 2 f3V 2 



2ir/3 



8tt 



(36) 



The counting field dependence is only through X\. This 
result obeys the FT since when A — > —A + t/3V, X\ is 
unchanged. Actually, already in Eqs. (fl7|) and (|26|) we 
can see the FT is satisfied since 



fi_A-M/3v( w ) = ^a(w) 



(37) 



This CGF possesses continuous lines of singularities 
for X\ <G [1, oo). As discussed in Sec. lU in the u— plane 
these become branch cuts (— oo,u_] and [u + ,0), where 
the branch points arc at u± = c^ l3V / 2 {x±^/x 2 — 1), with 
x = cosh(/3V72)(l - 1/T) - 1/T [see Figs. Q] (a) and 
(b)]. The zeroth-order CGF (f3"6")l takes a particular sim- 
ple form for a symmetric bridge Tl = at resonance 



eo = 0, for which the transmission is perfect T = 1 [see 
Eqs. J2ZI), QS), and (2)] 



J- (A) = zA(y + zA//?)/(27r) , 



(38) 



which describes Gaussian thermal fluctuations. On the 
other hand, when the bridge is extremely askew, Tr <C 
Tl or [eo | 3> T and thus T< 1 one may expand the CGF 
to obtain 



TV 



(- 



- 1 



- 1 



(39) 



2tt \l-e-P v eP v -l 
which is the bi-directional Poisson form. 

B. The phonon-induced part 



The noncquilibrium Luttingcr-Ward functional $ in 
Eq. (|18l) results from the coupling of the charge carriers 
to the vibrational modes, and therefore it is a functional 
of the Keldysh phonon Green function, 



d x (u) 



whose four components are 
2u> 



(40) 



7±± 



M =Rc 



(w + i0+) 2 - u; 2 
— z7rcoth — [5(o; — cjo) — <5(w + oJo)] 7 
d T± (io) = - 2m[5{u - uj a ) - 5(lu + w )] ?^ T H . (41) 

Here + is a positive infinitesimal, and n is the Bose 
distribution 



n ± (w) = ±- 



(42) 



o-o 




FIG. 4: Second-order diagrams, the Hartree (a) and Fock 
(b) terms. Diagram (c) is one of the fourth-order ones, and 
diagram (d) is one of sixth-order diagrams (a ring diagram). 
The full lines denote electron propagators, and the dashed 
ones phonon propagators. 

The functional $ can be expanded diagrammatically in 
powers of the electron-phonon coupling [sec Appendix [A] 
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for more details], or more precisely, in the small param- 
eter 



2^ 

Trr 2 



(43) 



Figure 0] shows the second-order diagrams [(a) and (b)] 
and two of higher-order diagrams [(c) and (d)]. Dia- 
gram (a) represents the Hartree term, Eq. (|A15|) . which 
is ignored below since it is independent of the phonon 
distribution. Diagram (b) represents the Fock term, 
Eq. (|A20|) . which depends on the Bose distribution. The 
actual noncquilibrium phonon distribution function can 
be obtained only by summing over infinite-order (in 
the elcctron-phonon coupling) diagrams. Therefore we 



carry out a re-summation of particular diagrams in the 
random-phase approximation (RPA) [see the discussion 
in Appendix Hand Eq. (|A21|) ] 



(44) 



* RPA (.9a) = -^J dukidet Dxiuj)- 1 , 
and then calculate up to second-order accuracy 



RPA/ 



(45) 



Here our expectation is that RPA would provide accurate 
result up to second order in 7 [see Eq. (fT8|)] except for 
the Hartree term, which we ignore. The dressed phonon 
Green function D\ is 



d(u) 



r 3 n A (w)T 3 = 



{u 2 - wg)/2w ■ 



M n+-( w ) 

-(^ 2 - Wo 2 )/2 Wo -n--H 



(46) 



with r 3 being the third Pauli matrix 



T3 



1 

-1 



(47) 



In Eq. (|46jl appears the Keldysh particle-hole propagator 
n, whose lesser/greater components are given by 



OI^M = £j dLj'g^{w' + U /2)g^{ U ' w/2) 



±t(A,-V) 



(48) 



r,r'=L,R 



where 



f_ f d l g^(u' + u/2)g^{u'-ul2) 
2tt J ^ n x (uj' + uj/2)n x (uj' -u/2) ' 

(49) 



(See Appendix [B] for more details.) The particle- hole 
propagator satisfies the relations 

n+-H = n^+(-w) , 

n+-( w ) = n-+(-^) , (50) 

with n^ _ (w) = 5Z rr , n+7^(u;). The analytic expression 
of Eq. (|49p in the extended wide-band limit is lengthy 
and is given by Eq. (|B3j) . For A = it attains a simple 
form 



with 



(51) 



(52) 



Relegating detailed derivations to Appendix |B] the 
causal and anti-causal components obey 



n++H + il--(cj) = &+-(«) + n-+( w ) + 0f ( W ) , 



n++M - nr-(w) = 2R e n fl (u;) + ^( w ) 



(53) 
(54) 



The definition of (f>f and its analytic expression in the 



extended wide-band limit are given in Eq. (|B7j) an d 
Eqs. (|B10|) and (|B11|) . respectively. One finds that 4>f is 
even in u. It turns out that <j) x may be safely neglected 
[see Eq. (|B15jl ]. 

Finally, the retarded component is given by exploiting 
the Kramcrs-Kronig relation 

n» = JL /^ n Q- + ^7 n ^). (55) 

' 2-k ] lj - lj' + i0+ y ' 

The particle-hole propagator is proportional to the small 
parameter g, Eq. (|43|) . 

Having determined the dressed phonon Green function 
Eq. (|4"B"]) . we now use it to obtain $ RPA , Eq. ([iijl. It is 
expedient to calculate its derivative with respect to A. By 
presenting the determinant of the dressed phonon Green 
function in the form [using Eqs. (|53[) and (|54|) ] 



dct D^fo) 



2w 



^-Ren J? (w) + .A A H 



(56) 



where 



^H=n A -+( w )n+-( w ) 

-i[n-+H + n+-H + 0fH] 2 , (57) 
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we obtain 

9$ RPA (gA) 
dX 



did 



dxA x (u) 



(58) 

The integral here is carried out to leading order in g. The 
integrand has four poles, which within that accuracy are 
given by 



±[wq + ReU R (tu )} ± iy/Ax(w ).. 



(59) 



where we have used A\(ui) = A\(— u). 

The locations of poles of the integrand in Eq. ([58| are 
complex functions of A. Consider them first in the limit 
A -> 0. Since 



A (uj ) 



I 2 ^ ga r r'Pl(u - Mr + Mr') 
\ r,r'=L,R J 



(2gpluj f 



(60) 



is a positive real number, then at small enough A, the 
real part of the squire root satisfies ReyAjJtJo) > and 
thus the integration yields 



up to 0{g). It follows that 



d\y/A\(u) ) 



(61) 



$< 2 > (A) = y/Ax(uj ) - 2gf&o , (62) 

which satisfies the normalization condition <I>( 2 )(0) = 0. 
Away from the origin A = 0, Tie Ax (u>o) may be nega- 
tive. When this is the case we obtain 



l> (2) (A) = - y/Ax(u ) - 2gp 2 u J o ■ 



(63) 



Hence a branch cut appears at A satisfying 
Rc^aM = 0. 

In order to carry out a large-deviation analysis it suf- 
fices to consider the imaginary axis of the A—plane, since 
there the rate function ((SJ is real. On the imaginary axis, 
A\(u)o) is a rca l function. When A\(u>o) > the imagi- 
nary part of <I>( 2 ' is zero and thus the imaginary axis can 
serve as the steepest patb4& of the integral Eq. Q. When 
A\(uio) < 0, there appears a branch cut on the imaginary 
axis. One might have thought that the branch cut would 
be detrimental to the saddle-point approximation of the 
integral in Eq. (|4")l. Here we point out that a complex in- 
tegration along a path encircling the branch cut oscillates 
rapidly in the r — ¥ oo limit and thus would be averaged 
out. Therefore we will neglect the contribution from the 
branch cut. 

Equations (|62|) and ([63]) are the main results of this 
paper. Similar expressions have been reached by Haupt 
et. al. in Refs.[l9t these authors concentrated on the first 
few cumulants. Here we obtain analytical expressions for 



the CGF; those enable us to examine its singularities and 
to fully analyze the rate function itself, as is detailed in 
the next Section. 

Finally, we verify that our results obey the FT. Ex- 
ploiting the detailed-balance relation, 



a /9(w-jUr+/V) 



(64) 



(65) 
(66) 



n-t_A+w( w ) = n rrt aH = n+7»e 

and the relations 

n:+ +w M = n+-(c)e^ ; 

we find 

A-\+ipv(u) = A\(w) . 
Therefore Eqs. flfE) and ([63]) satisfy the FT, 

¥ 2 \-X + i(5V) = ¥ 2 \X). (67) 

IV. RESULTS AND DISCUSSION 

In the following we confine ourselves to a symmetric 
junction, a = 1 (IY = Tr), and thus the normalized 
density of states on the dot at the Fermi level dominates 
the transmission probability, T — po [see Eqs. (|24|) . (|27j). 
(f2"8")h and ([3T|) ]. The numerical results are all obtained for 
the electron-phonon coupling constant, Eq. (|43"]) , g = 0.1, 
unless otherwise specified. 

A. Average current and noise 

Figure [5] (a) shows the source-drain bias voltage depen- 
dence of the current ((/)) at a finite temperature /3ujq = 10 
(solid lines) and at zero temperature (dashed lines) for a 
perfect, T = 1, and a relatively weak, T = 0.5, transmis- 
sion probabilities. At perfect transmission the current is 
suppressed above the threshold, \V\ > wq, because elec- 
trons are inelastically backscattered by phonons. When 
the transmission is weak, T = 0.5, the current is slightly 
enhanced above the threshold. These results are consis- 
tent with previous ones^ A finite temperature tends to 
smear the kink structure of the T = 1 curve; it affects far 
less the average current at weak transmission, T = 0.5, 
where the solid and dashed lines almost overlap. 

Figure [5] (b) depicts the current noise ((I 2 )). At per- 
fect transmission the noise is absent below the thresh- 
old \V\ < ljo at zero temperature. Thermal fluctua- 
tions which arise at finite temperatures induce additional 
noise below the threshold. Although the current is sup- 
pressed above the threshold |y| = ojq by the inelastic 
phonon scattering, the noise is significantly enhanced. 
This indicates that inelastic phonon scattering broadens 
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the probability distribution of the current. In the case of 
a weak transmission, T = 0.5, the temperature effect is 
less dramatic-the noise is simply enhanced. 



3 
P4 
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FIG. 5: The source-drain bias voltage dependence of the 
current (a) and the current noise (b) for T = 1 and 0.5. 
Solid lmes-/3u}o = 10; dashed lines-zero temperature. The 
vertical axes are normalized by loq/Rk, where 71k = 2n is the 
resistance quantum. 



B. Singularities and the rate function 

At zero temperature, it is possible to obtain a simpler 
form for the scaled CGF, which is useful for finding its 
singularities. For positive voltages V > 0, the explicit 
form of the electronic part of the CGF, Eq. (|3"6"|) . reads 

FoW = 77- ln "> u = l + T(u-l), (68) 

where u is defined in Eq. ([3]). The function A\ which 
determines the phonon part is given by 

Ax(u>o) -4 5 2 T 2 [(1 - l/u)(T - l)V + T^o] 2 , 

for < V < uo and 

AxM =g 2 T 2 {[V(l -T + u 2 (2T - 1)) + (T- 1 

+ fi(2 - u + 2T(u - 1))) w ] 2 - uT 2 (V 



(69) 



- lu )[2ulu + u(V + (2« - 1) wo)]}/« , (70) 

for V > Ldo. Below, we investigate the analytic properties 
of the phonon-induced part of the CGF. 



1. Elastic phonon scattering 

As at zero temperature phonons cannot be excited 
when < V < loq, and electron transport at such volt- 



ages is hence affected only by elastic phonon scatter- 
ing. Figure |5] (a) shows schematically the square-root 
branch cut of # (2) , Eqs. (JUJ) and {B3J). We find that 
the u — plane is separated into two by the the brach cut. 
This branch cut intersects the real axis at Uo = 1 — 1 jT 
[uq is indicated by empty dots in Fig. [5] (a)] and u\ = 
(T - 1)(V + wo)/[V(T-l)+Tw ]. Around these points, 
A\ can be expanded as 



A> 



(u - u ) 2 



Ax w g 



AT\Tuj q + (T -l)V) 4 
V 2 (T-1) 2 



(u - ui) z 



(71) 



(72) 



Upon sweeping the transmission T from to 1, uq in- 
creases from — oo to 0, while u\ increases but at 7c = 
V/(V + ujq) jumps from +oo to — oo. Accordingly, we 
may define two regimes, or phases 2 ^, I (T < 7c) an d II 
(T > 7c) as indicated in Fig. [6] (a). This classification 
roughly captures the behavior of the average current and 
the current noise. Figures[6](b) and (c) depict the correc- 
tions induced by the electron-phonon interaction in the 
average current and in the current noise, respectively, 



<(n> P h = 



9™$( 2 ) 



x=o 



(73) 



We find that the electron-phonon interaction always in- 
creases the average current under the conditions adopted 
here [Fig. |6] (b)]. On the other hand, the noise can be ei- 
ther enhanced or suppressed, depending on which regime 
the transmission is in [Fig. [B] (c)]. 

In regime I, one of the intersection points, Ui, is on the 
positive real axis, outside the unit circle \u\ = 1 [the left 
panel of Fig. [5] (a)]. Therefore, in the A— plane, there 
is a nonanalytic point on the positive iX— axis, which 
induces a weak non-convexity of the CGF as shown in 
Fig. [7] (a) [The non-convex region is indicated by an ar- 
row in Fig. [7] (a). For comparison, we also plot the g = 
case for which the CGF is convex, shown by the dot- 
ted line.]. Figure [7] (b) exhibits the rate function, the 
Legendre transform of the CGF, Eq. ([8]). The rate func- 
tion is multi- valued around I/((I))o ~ 3 (the thick line 
in the figure) because of the non-convexity of the CGF. 
As seen from its definition Eq. ((HJ, the rate function is 
obtained by choosing the A that maximizes the expres- 
sion in the curly brackets there. Then, similar to the way 
a first-order phase transition manifests itself in thermo- 
dynamics, a kink appears in the rate function. We note 
that the location of the kink does not coincide with that 
of the peak; the peak of the rate function is at A = 
where 1 = and / = ((/)) . The physical consequence is 
that the elastic phonon scattering broadens the probabil- 
ity distribution by enhancing the probability to observe 
currents larger than the average value. The kink is a con- 
sequence of the non-convexity of the CGF on the positive 
real axis of the u— plane, a feature which is absent in the 
nonintcracting electron case, and is a signature of the 
electron-phonon correlation. 
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In regime II the two intersection points are on the neg- 
ative real axis of the u— plane [the right panel of Fig. [5] 
(a)] . In this case, the CGF can be reduced to that of non- 
interacting electron transport^ The CGF and the rate 
function pertaining to this case are plotted in Figs. [7] (c) 
and (d) . The CGF is convex and the corresponding rate 
function is concave. In this regime the peak position is 
shifted from that for g = [dashed line] . The elastic scat- 
tering by the phonons can either broaden or shrink the 
width of the rate function depending on the transmission 
probability as we deduce from Fig. [5] (c). 
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FIG. 6: (a) Schematic picture of the branch cut of <J>' 2 - ) in the 
complex it— plane below threshold, < V < ujq- The dashed 
circles are the unit ones. The corrections induced by elas- 
tic electron-phonon scattering to the average current and the 
current noise are portrayed in panel (b) and (c), respectively. 
The bias voltage is V/uJo — 0.5. See text for the significance 
of To. 
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FIG. 7: The scaled CGF (a) and the rate function (b) in 
regime I (7~ = 0.1). Solid lines are for g — 0.1 and dashed lines 
are for g = 0. The bias voltage is below threshold, V/ljq = 0.5. 
The CGF is plotted as a function of the counting field on the 
imaginary axis. The axes are normalized by {{/))o = TV/Rk- 
The scaled CGF (c) and the rate function (d) in regime II 
(T = 0.5). 



2. Inelastic phonon scattering filled dots in Fig. [5] (a)]. In addition there is another 

point, M = 1 — 1 /T [the empty dots in Fig. [5] (a)] , around 

Above threshold V > uq, the inelastic phonon scat- which A * can be ex P an ded 

tering becomes possible. The analytic properties of the 2V(T — l)(wo — V) 

scaled CGF in this regime depend on the bias voltage. A\~g — - . (74) 

The branch cuts in the u— plane are schematically shown 

in Fig. [5] (a). Three branch points, Ux, u 2 , and u 3 can The branch point u\ is always on the negative real axis 

be obtained by searching for the roots of A\(u) = [the such that u\ < Uq. From the positions of the other two 
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branch points, 1*2 and U3, we identify three regimes, see 
Fig. [5] (a). In regime I the two branch points are on the 
positive real axis outside the unit circle [dashed circles in 
Fig. |S] (a) are unit circles]. As T increases u-i and 1*3 ap- 
proach one another and meet at T = T- ■ Then in regime 
II, the two branch points are located symmetrically off 
the real axis. Upon further increasing T, and 113 move 
in the complex plane and at T = 7+ they meet on the 
real axis again. In regime III, the two branch points are 
on the positive real axis inside the unit circle. 

In Figs. [5] (a) and (b) we plot the CGF and the rate 
function pertaining to regime I. The overall tendencies 
are similar to those found in regime I below threshold, 
< V < l)q. In the shaded area of Fig. (a) the CGF is 
nonanalytic and non-convex. This gives rise to a stronger 
kink structure in the rate function [Fig. [9] (b)]. Figures 
M (c) and (d) present the CGF and the rate function in 
regime II. The overall tendency is again similar to those 
of regime II below threshold. However, because of the 
non-analyticities off the real axis, the statistics would 
not be reduced to that of noninteracting electrons. 

We note that our identification of the three regimes 
roughly captures the behavior of the corrections to the 
current and the noise induced by the electron-phonon in- 
teraction [Figs. [H (b) and (c)], which oscillate as a func- 
tion of the transmission probability^^ 
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The analysis of the CGF and the rate function in 
regime III [see Fig. [5] is rather subtle. At zero temper- 
ature, there develops a non-convex region in the CGF, 
and the origin A = might enter it. When this happens, 
second-order perturbation fails since the rate function 
does not satisfy the relation 1(1 = ((I)) ) = [this relation 
is related to the normalization condition J-(X = 0) = 0]. 
We therefore study regime III at finite temperatures, tak- 
ing as an example (3wo = 10, and show that the FT is 
crucial for obtaining a physically reasonable result. Note 
that when the symmetry holds, the FT ^ is also 
preserved within the large-deviation analysis^ 



FIG. 8: (a) Schematic picture of the branch cut of $ ( - 2 ^ in the 
complex u— plane above threshold V > cjo. The dashed circles 
are the unit ones, (b) The corrections induced by the electron- 
phonon coupling to the current and (c) to the current noise. 
The bias voltage is V/ujq = 1.5. The boundaries between I and 
II and between II and III are at T- ~ 0.551 and 7+ ~ 0.908, 
respectively. Regime II roughly corresponds to the negative- 
noise correction region. 



1(1) = max {£1 - + ifiV)} 

= max{-r/-/3/^-^Hn} 
=1(-I) - P IV . 



(75) 



Figure [TU] exhibits the CGF and the rate function at 
perfect transmission. For comparison, we plot the corre- 
sponding curves for noninteracting electrons, Eq. (|38j) 
[dashed lines]. As we have already noted when dis- 
cussing the current noise, Fig. [5] (b), the width of the 
rate function is enhanced by inelastic phonon scattering 
[Fig. M (b)]. The CGF obeys the FT, Eq. ©, and 
the curves are symmetric around the dot-dashed verti- 
cal line at iX = -0V/2 [Fig. [TO] (a)]. The peak of the 



probability distribution is shifted in the negative direc- 
tion and the probability to find large current fluctuations 
is suppressed as compared with the noninteracting case 
[Fig. [TO] (b)]. In the shaded area of Fig. [TO] (a), the CGF 
is non-analytic and non-convex. Correspondingly, the 
rate function has a non-differcntiablc point at / = 0, see 
Fig. [TO] (b). As a result, although the probability to ob- 
serve currents smaller than the average value is enhanced 
by the inelastic phonon scattering, the probability to find 
negative currents I < is strongly suppressed. This is 
consistent with the FT ([6]), which states that although 
thermal agitations generate current flowing in the op- 
posite direction to the source-drain bias, that current is 
exponentially suppressed at low temperatures. 
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FIG. 10: (a) The CGF and (b) the rate function for T = 1, 
V/ujo = 1-5, and /3cjo = 10. The dotted (solid) lines show 
results with (g — 0) and without (g = 0.1) electron-phonon 
interaction, respectively. In the shaded area of panel (a) the 
CGF for interacting case is non-analytic and non-convex, re- 
sulting in a non-differentiable point of the rate function at 
1 = (b). Axes are normalized by ((/))o = V/Rk- 
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FIG. 9: The CGF (a) and the rate function (b) in I (T = 0.1). 
The solid lines are for g = 0.1 and the dashed ones are for g = 
0. The bias voltage is above threshold, V/uiq = 1.5. Panels 
(c) and (d) show the CGF and the rate function, respectively, 
pertaining to II (T = 0.6). Axes notations are the same as in 

Fig. m 



D. Discussion 

Recently, Kumar et al — suggested to identify three 
regimes, based on the sign of the correction of the noise 
induced by the coupling to the phonons, as depicted in 
Fig. [5] (c) . Our classification roughly agrees with that of 
Ref. lf)|, save for the transition points. The ones quoted 



there, T£ umal = 1/2 ± 1/(2^2) « 0.85,0.15, arc dif- 
ferent from ours and moreover, our transition points 7± 
depend on the bias voltage. Hence, there is no one-to-one 
correspondence between the classification of Ref. and 
ours. Since oscillations in higher cumulants are ubiqui- 
tou a 20 ' 53 ' 55 and dominated by locations of singularities, 
we think it would be legitimate to utilize the location 
distribution of singularities for the classification^—. 

A comment on the validity of second-order pertur- 
bation theory is called for. Flindt et. al^ analyzed 
the FCS of sequential transport through a quantum dot 
containing two levels^ They assume identical incoming 
rates but different outgoing rates for the two levels and 
analyzed the approximate CGF, obtained by expanding 
in the ratio of the two outgoing rates, taken to be a small 
parameter. They noticed a peculiar behavior: The agree- 
ment between the cumulants obtained by differentiating 
the approximate CGF (with respect to the counting field) 
and those derived by differentiating the CGF before ex- 
panding it was good for the low-order cumulants, but 
failed completely for the higher-order ones. In our case, 
the CGF is obtained as an expansion in the electron- 
phonon coupling and thus it is plausible that the unphys- 
ical results which we encountered around regime III for 
V > wo at zero temperature, may have a similar origin 
to the apparent discrepancy reported in Ref. l56l . This 
fault may be resolved by accounting for all orders in the 
electron-phonon coupling. However, getting analytical 
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results seems to be technically complicated and almost 
inevitably requires numerical methods^ 

Another example of the FCS for interacting electron 
systems is found in transport through a quantum dot 
in the Kondo regime. Recently, Sakano et. al. have 
calculated the CGF for the SU(N)— impurity Anderson 
model, — using the renormalized perturbation theory. 
Their result is exact up to the cubic order in the source- 
drain bias voltage 0(V 3 ) for a particle- hole symmetric 
case, and is quadratic in u. This is interpreted as the 
sum of the CGF for single-particle transfer and that for 
two-particle transfer. In our case, the phonon part of the 
CGF at zero temperature, to the leading order in V, can 
be obtained from Eqs. (|62|) and (|69|) in the form 



i>(2) 



n=i v ' 

2gT(T-l)(u-l) 



U — Uq 



(76) 



where uq = 1 — 1/T. One may interpret the first right- 
hand side of Eq. ()76|) as a sum of independent n-electron 
transfers. However, the second line indicates that the 
CGF is non-analytic on the negative real axis of the 
u — plane at Uq < 0. Therefore, the electron transfer 
statistics can be reduced to that of noninteracting elec- 
trons^ This example suggests subtleties in interpreting 
the CGF for interacting electrons. 



V. SUMMARY 

We have investigated the full-counting statistics of cur- 
rents mediated by clastic and inelastic elcctron-phonon 
scattering. In the extended wide-band limit, we obtained 
analytical expressions for the cumulant-generating func- 
tion, accurate up to the second order in the electron- 
phonon coupling. Our analytic expressions are applica- 
ble for finite temperatures and bias voltages and satisfy 
the fluctuation theorem. Using those we analyzed the 
locations of singularities of the CGF. The singularities 
are symmetrically distributed in A— plane because of the 
fluctuation theorem. 

The singularities in u— plane, which appear because of 
the electron-phonon interaction, classify specific regimes 
in which the dependence of the electron transfer statis- 
tics, the current and the current noise, on the bare trans- 
mission is distinct. For small transmission probabilities 
we find singularities of the CGF on the positive real axis 
satisfying u > 1. Around the singularities, the CGF is 
non-convex, which results in a kink of the rate function. 
Such a kink, derived within the large-deviation analysis 
resembles a first-order phase transition in thermodynam- 
ics. It signifies the tendency of the phonon scattering to 
enhance the probability to find currents larger than the 
average value. 

When the bias voltage is larger than the phonon fre- 
quency, V > u>o, we find singularities in < u < 1 



around perfect transmission. This results in a kink at 
I = of the rate function. Consequently the rate func- 
tion is broadened because the probability to find currents 
smaller than the average value is enhanced. In the inter- 
mediate transmission probability region a negative cor- 
rection in the current noise is induced by the inelastic 
electron-phonon scattering. 

The kink structures in the rate function character- 
izes the electron-phonon interactions. Although mea- 
surements of the rate function would be technically 
demanding, the FCS can be in principle obtained 
experimentally^ giving hope that our predictions could 
be put to test. 
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Appendix A: The Luttinger-Ward potential 

The way to construct a partition function based 
on the self energy is to exploit the Luttinger-Ward 
functional approach^ or the self-consistent ^-derivable 
approximation.— This method can be straightforwardly 
extended to a nonequilibrium situational The saddle- 
point approximation for the CGF can also be constructed 
by this approach^ The underlying idea is the observa- 
tion that by introducing the Luttinger-Ward functional 

which includes all skeleton diagrams, the total gener- 
ating functional can be formally written as 

J r (X)= lim -fTr[lnG- 1 ]+Tr[S(G)G]-$(G) N ) , 

T— >oo T \ J 

(Al) 

where the trace and the product mean integrations over 
time along the Keldysh contour. The Green function G, 
also defined on the Keldysh contour, is 



G- 1 (1,2)=. 9 - 1 (1,2)-S(1,2;G) , 



(A2) 



where the arguments 1, 2, • • • stand for t\, £2, • • • • Here 
gx is the Green function of the noninteracting electrons 



ffA (l,2) = -i(T K c (l) I 4(2) I ) 



(A3) 



given explicitly in Eqs. (f2"2"j) and (|23p . The self energy is 
a functional of G as well as of the phonon Green function 
d, 

d(l, 2) = -i (T K (6(l)i + 6 t (l)i)(6(2)i + 6 t (2) I )) , (A4) 
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given in Eqs. (|40|) and (|4Tj) . The functional derivative 
of the nonequilibrium Luttinger-Ward functional $ gives 
the self energy, 



£(1,2; G) 



<5G(2,1) 



(A5) 



Both functions G and £ depend implicitly on the count- 
ing field only through g\. 

By differentiating Eq. (|A1[) with respect to the count- 
ing field A and using Eq. (|A5[) one obtains 



<M _ _ to I Tr f G %l 

CtA T^oo T \ C/A 



where we have used the relation 

d$ _ / <5$ <9G 
rfA ~~ Tr l 5G "aX 



(A6) 



(A7) 



It should be emphasized that in this formulation the self 
energy has to be determined self-consistently in order 
to satisfy conservation laws4£ To second order in the 
electron-phonon coupling 7, Eq. (|A6[) becomes 



dA 



lim -Tr [ g x - 

T->oo T \ OX 



lim -Tr ( gx^(gx)gx^- 

T->00 T V CM 



where we have replaced S(G) by Y,(g\), since it is already 
0( 7 2 ). Upon using the identity 



dgx 
dX 



-gx- 



dgl 1 

dX 



-fix 



(A9) 



and exploiting Eq. (|A5]) [in its 0(-y 2 ) form] Eq. (|A8|) 
becomes 

— ps hm -— Tr ln.g A ] - hm — — , (AlO) 

d\ t^oo t dX t->oo r dA 

and consequently 

T(X) = lim -(Trpn^ 1 ]) - & 2 \X) . (All) 



Comparing this expression with the original one, Eq. 
(|Aip . one notes that the term depending explicitly on 
the self energy has disappeared. The lowest-order scaled 
Luttinger-ward potential $( 2 )(A) is obtained by expand- 
ing $(gx) (which depends on gx and not on G) up to 
C( 7 2 ) [see Eq. (18])]. 

For brevity, the calculation above is presented in the 
time domain. One may also switch to the frequency rep- 
resentation in which Eq. (|A6[) reads 



AF(A) 
dX 



1 

2^ 



ox 
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The Pauli matrix T3, Eq. (j47J), appears once we project 
the time from the Keldysh contour on the real time axis, 



,t/2 ,t/2 

dt= dt + - dt. 

J-t/2 J-t/2 



(A12) 



where t± e K± (see Fig. [3]). This leads to Eqs. flU) and 
(jTTJ) in the main text . 

It remains to calculate &(g\) up to 0( 7 2 ). This is 
accomplished by expanding it perturbatively in 7. The 
zeroth-ordcr term is just a constant, 



* (0) (5a) 



lim — Trln[d] 



(A13) 



independent of the counting field. The diagrams con- 
stituting the second order are depicted in Fig. [4] The 
Hartree term [Fig. @] (a)] is 

* H (.9a) = -^ Um - / dld2g x (l,l)d(l,2)g x (2,2) 

ss'd ss '(0) Jduigi'iuii J du 2 gi s \u 2 ). 

(A14) 



S,S =+,— 



(A8) Inserting Eqs. (J33J) and flU) into Eq. (|AT4)) yields 



^0 \j 



n A (w) ( w -e ) 2 + r 2 /4 



Oa(w) 



, (A15) 



where the small parameter p is given in Eq. (|43j) . Note 
the appearance of the distribution Eq. (f2"Tj) in the form 
/ r + (w) — 1/2, resulting from the definition of the step 
function as 0(0) = 1/2 in the continuous notation.— By 
using Eq. (|3"T|) . one can check that the FT is fulfilled, 
i> H (g_A+i/3v) = i> H (<7 A ). In the main text, we neglect 
the Hartree term since it does not depend on the phonon 
distribution. 

The Fock term [Fig. H (b)] reads 



® F (gx)=^ lim - / dld2«ft(l,2)d(l,2) 5A (2,l) 

= -- lim - / dld2d(l,2)II A (2,l) , (A16) 

where we have introduced the particle-hole propagator, 

n, 



n A (l,2) = -z 7 V(l,2) 5 A(2,l) 
Adopting the form (|A16[) . we find 



(A17) 
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* F (3a) 



1 f du) 
2^ 



d ++ (uj)U+ + (uj) +d~~(w)nr _ (a;) - <i" + (w)n+"(w) - <i + -(w)II^ + (tj) 



(A18) 



Here a constant should be added to keep the normalization condition $ F (<7o) = 0. Inserting Eqs. (|53|) and ((54)) . we 
obtain 



* F (5a) 



1 [ duj 



d~ + (oj) 



n+-(w) + n+-H 



n+"(o;) +d+-( w ) 



2 7 2tt 

o 0A M + o H 



(A19) 



Using relations (f5T)|) and the corresponding ones for the components of the phonon Green function, e.g., d + (w) = 
h (— w), yields 



ifl x + (suj ) + ill + {soj ) 



iH x + (sujq ) ) + - coth 



2 



^ (wo) 



duj cjo <t>\ (w) 

27T CJ 2 — CJq 



(A20) 



where P means Cauchy principle value. By using 
Eqs. ([65j and (|B13|) . one may verify that the FT is 
obeyed, $ F (g_ A+t ,gy) = $ F (ff A )- 

The Fock term depends on the equilibrium phonon dis- 
tribution n + ) 19 ' 21 which suggests that a rcsummation of 
an infinite series of diagrams is needed in order to ac- 
count for the noncquilibrium phonon distribution. We 
carry out this summation within the random-phase ap- 
proximation, by summing over all ring diagrams [see di- 
agrams (c) and (d) in Fig. [4]. The RPA is also known 
to be relevant for the AC conductance^ In this way, we 
obtain the functional 



' PA (3A) = i lim -Trln[l-dn A ] + i>(°)( 5A ) 

Z t — f oo 7- 



1 



1 



— lim — TrlnD 

2 t— >oo T 



A ' 



(A21) 



where we have also included the zeroth-order term, Eq. 
(IA13[) and the full phonon propagator, D\, 



d: 



d- 1 -m 



(A22) 



Equation (|A21[) yields Eq. (|44p in Fourier space. 



This RPA can be also formulated using the Keldysh 
path-integral approach within the saddle-point approx- 
imation^ and further accounting for the Gaussian- 
fluctuation correction around it. 



Appendix B: Analytic expressions for the 
particle-hole propagator 

Here we obtain analytic expressions for the four com- 
ponents of the particle-hole propagator, Eqs. (|49|) . in the 
extended wide-band limit J^. Since the lesser and greater 
components are related to one another [see Eqs. ([5U| ] it 
suffices to compute the greater function, which reads 

in rr J x {uj) ss ga rr >p / du , ■ , r 7 _ , , (Bl) 



Oa(w+)Oa(w-) 
where the small parameter g is given in Eq. (|43[) . and 

lo ± = w'±Ui/2 , (B2) 
A lengthy calculation yields 



ffl-+ A H = ^^e^" +s ^)/ 2 ^(z 4 + e-^-^/ 2 )(z 4 +e^^'^/ 2 )logz i n^— . (B3) 



where s r = ±1 for r = L/R, 

zi — e^ Z A+ , z 2 — e^ Z A _ , z 3 — e z A+ , z 4 — c z A _ , 



r 



and Z" A ± are given by Eqs. (|3"3")l . The calculation of tionsn ++ ±n . The only A— dependence of the casual 
n ++ and n is facilitated by considering the combina- 
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and anti-casual electronic Green functions is contained 
in their denominator, il\ [see Eqs. (f2U)) ]. Therefore we 
may write 



Upon using the relations 



7 2 [^.g$ ± (u, + )g^(uj-) 



±±t 



(B4) 



.9o ++ = ,9o + + 9 R = % + + 9 A , 
9o~ = 9o + - 9 R = 9o~ - 9 A . 



(B5) 



where g R ' A are the retarded and advanced Green func- 
tions, 



1 



(B6) 



uj — e ± iT/2 ' 

we find Eq. ([S3J with 

2 r I 
ifli>f (w) =— / do/( ; : ; - ] 

x [g R (u + )g R (uj_) + g A (u J+ )g A (u J _)} , (B7) 
where the relation 



du/g R l A (u.)g R / A (w_)=Q 



(B8) 



has been used. Since fl\ = 1 for \ui\ ^ max(|U|, 1//3) 
the integral in Eq. (|B7|) is bounded, and therefore in 
the extended wide-band limits the terms in the square 
brackets there can be replaced by 



2Rc 5 i? (0) 2 



yielding 



i<t> s x (u) 



7 [i-2/7 ]/r 2 (|T/|,fc B r, Wo <r) 

2/e 2 (|7|,fer > w ,r«|eo|) 

(B9) 



2 ffPo [l-2po]^ (|V|,fcBT,w «r) 

(|V|,fcBT,wo,r«|co|) ' 

(BIO) 



where 
S= I doj' 



- 1 



n 



i \ s.s' — ± 



(BH) 



Turning now to the combination II ++ — II , wc use 
Eqs. (|B5P and the relation g K = g^ + + <?q for the 
Keldysh component of the Green function, to obtain 
Eq. ([54")) . Equation ([55| is derived from the relation 

n « n= y /^ g K K)g A (^-) + .g fl K).9 K (^-) 

[UJ> 2m J ^ 2 
In Eq. (|5"4"|) . we obtain, 



2 /* i 

= X. / d^'f — * — _ - 1) 



x Im(. 9 K K)/( W _) + g R (u; + )g K (u_)) . (B12) 
Using Eq. (J3TJ) , we can verify the symmetry 



4c x+iflv (u) = ^(«). 



(B13) 



One can now convince oneself that in the extended wide- 
band limit, in which Re[g R (oj)} « — e /[eo + T 2 /4] and 

5 K M = - r E f Potanh[/3( W - Mr )/2] (B14) 



is C(l/r), </>^ may be safely neglected, since 

^ K 1 l/(r|e |) (|U|,fc B T, Wo ,r« |e |) ' [mb) 
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